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Abstract 

We introduce a new model of random layered media, extending the Matheron- 
de Marsily model: Here we allow for the flows to change in time. For such layered 
structures, we solve exactly the equations of motion for single particles, and also for 
polymers modelled as Rouse chains. The results show a rich variety of dynamical 
patterns. 

PACS numbers: 36.20.-r, 47.55.-t, 05.60.-k 



1 Introduction 

The dynamics of particles and particle assemblies in random force fields is a subject of 
intense current interest (see Ref. [1-3] and references therein). Such fields often lead to 
qualitative deviations from simple dynamic patterns, and result in the anomalous diffusion 
of the particles involved. A standard prototype model is that of Matheron and de Marsily 
flows (MdM) [1-10], originally designed to describe the transport of a solute in porous 
media. Similar mathematical forms arise when treating electrons in random potentials or 
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spin depolarization in random fields [10]. In the MdM model the environment is viewed 
as consisting of layers of force fields [4] . The particles' mean square displacement (msd) 
in the direction parallel to the flows, say, along the Y-axis, < Y 2 (t) >, can be evaluated 
exactly; it has been shown that it grows as < Y 2 {t) >~ t 3 / 2 [4,5,8], t being the time. 
Recently the dynamics of Rouse polymer chains in time-independent MdM flows were 
also discussed and several exact results have been derived [1, 2]. In particular, it has been 
shown that the msd of some tagged bead of the chain may display different dynamical 
regimes depending on whether the time of observation is shorter or longer than the so- 
called Rouse time £r. Now t R goes as N 2 , where N is the degree of polymerisation, i.e. 
the chain's length [11-13]. For times shorter than t R a sub-ballistic law < Y 2 (i) >~ t 7 / 4 
was derived, while it was shown that for times longer than t R the t 3//2 -dependence is 
restored, the prefactor being a growing function of N. Thus a tagged monomer of a 
polymer immersed in MdM flows moves in general faster than an individual particle (i.e. 
a chain with N — 1), a fact which has its physical explanation [1,2]. 

These results concern, however, velocity fields whose directions and magnitudes are 
random in space but fixed in time: we speak of quenched disorder. For random diver- 
genceless velocity flows, for which the flows change with time, we show in the following 
that new results emerge. 

In the present paper we present exact dynamical solutions both for individual particles 
and for long Rouse chains moving in time-dependent, random layered media. Distinct 
from the previous works [1,2,4,5,8], we now allow the velocities in the layers to change 
randomly with time. We thus take into account the fact that the environment itself may 
be subject to dynamical processes, thereby introducing new time-scales into the problem. 
By this the time fluctuations of the flow give rise to interesting, novel dynamical behaviors. 
In particular, we show that at intermediate times (t <C t R ) the msd of any tagged bead of 
a polymer chain follows < Y 2 (t) >~ t 3 / 4 , i.e. the msd moves faster than the Rouse law 
< Y 2 {t) >~ t x l 2 which holds in the flow-free case [12, 13]. On the other hand, at times 
t 3> tR, we recover a diffusive behavior, in which however, the effective diffusion constant 
scales as iV -1 / 2 for sufficiently long chains, i.e. vanishes with N at a slower rate than the 
Rouse diffusion constant, which goes as 1/N. 

The paper is structured as follows: In Section II we formulate the model and introduce 
the basic notation. In Section III we analyse the dynamics of an individual particle subject 
to time-dependent MdM flows. Section IV is devoted to the analysis of the dynamics of 
a tagged bead of an infinitely long Rouse chain, while Section V deals with effects typical 
for finite chains. We conclude with a discussion and a brief summary of results in Section 
VI. 
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2 The model 



A standard polymer model due to Rouse [11] consists in viewing the macromolecule as 
a series of N beads, linearly connected by harmonic springs. Its dynamics in solution, 
but subject to external force fields and excluding hydrodynamic interactions and steric 
hindrances, is given by the Langevin equation [11-13]: 

= K(R n+1 + R n _! - 2R n ) + F(R n ,t) + n(n,t), (1) 

for 1 < n < N, complemented by R = Ri and R^ = Rjv+i- In this equation ( is 
the coefficient of friction, R n is the position of the nth monomer, K = ?>k^T jb 2 is the 
spring constant accounting for the Hookean interaction of the monomers, ksT being the 
temperature multiplied by the Boltzmann constant and b the average distance between 
the beads, and 77(71, t) represents the thermal noise due to interactions with the solvent. In 
Eq.(l) the F(R n , t) denote the extra forces due to the imposed velocity fields. Equation (1) 
is a simplified description of the dynamics of a polymer in ^-solutions [12, 13]. 

Treating the index n as continuous leads to the replacement of the discrete term 
^R„ + i + R„_i — 2R n j by the Laplacian operator d 2 H n /dn 2 and of R = R x and H N = 
Rjv+i by the Rouse boundary conditions [12]: 

<9R 

— - = 0, for n = and n = N. (2) 

on 

Next, we specify the forces entering the right- hand- side (rhs) of Eq.(l). We let, as 
usual, the thermal noise be Gaussian with zero mean; due to the fluctuation-dissipation 
theorem the second moment takes the form: 



i] a (n, t)np(n', t, ) — 2ksTC > 5 a ^5 n ^ n i S(t-t'), (3) 

In Eq.(3) we let the Greek indices denote the cartesian components. The Matheron-de 
Marsily (MdM) model is obtained by assuming layered flows, say by taking the forces to 
be along the F-axis, but to depend solely on the X-component, see e.g. Fig. (1) in [1]: 

F(R n ,t) = (0,/(X n ,t),0). (4) 

This choice of F models a layered medium along the X-axis; each layer has a random (but 
for all points of the layer fixed) velocity pointing along the Y-axis. Due to the uncoupling 
of the different components in the Rouse-model, the dynamics in the transversal plane, 
(i.e. the X- and Z-components) are not influenced by the presence of the velocity fields. 
However, as we shall see, the motion parallel to the F-axis is dramatically changed. 
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The random force F(R n ,t), (i.e. its F-component) is assumed to be Gaussian and 
zero-centered. For the second moment we take 

(f(X, t)f(X>, f )> = ^S(X - xy-^ (5) 

Note that the strength of the random velocity field depends on the constants T and A. 
Equation (5) defines a changing environment with a short term memory, exponentially 
decreasing on the time scale T; we will call T the renewal time. In the limit r — > oo and 
keeping the ratio A/T fixed we recover the standard quenched MdM model. We shall 
denote this limiting procedure as the quenched limit in the rest of the paper. On the 
other hand, when r — > 0, exp(— |t|/r)/(2r) tends to the Dirac 8(t) distribution; this limit 
allows us to study an environment without memory. We distinguish averages over thermal 
histories from averages over configurations of the velocity field by using an overbar for 
the former and angular brackets for the latter. 



3 Dynamics of a single particle 

To fix the ideas, we start with the case of a single bead in time-dependent MdM flows. 
Then N — 1 and K — 0. Because of the decoupling of the different coordinates, the X- 
and Z-components of the particle's displacement obey the force-free Langevin equation 

. dX , . , . 

C^ = ^(t), (6) 

with T) x (t) being the X-component of the thermal noise. Equation (6) yields for the mean 
square displacement (msd) averaged over thermal histories the usual diffusion result: 

knT 



X\t) = 2-^- = 2Dt, (7) 

Through the last relation we have introduced the "bare" diffusion constant D. For the 
F-component we have to solve the equation 

C^=Vy(t) + f(X(t),t) (8) 

with f(X(t),t) being the random function of Eq.(5). The solution of Eq.(8) is: 

Y(t) = y*dT( Vy (T)+f(X(T),T)), (9) 
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where we set Y(0) = 0. Being interested in the mean square displacement we have: 

=2Dt+^ f dh fdhifiXitJ^fiXih)^)), (10) 
s Jo Jo 

by noting that the different sources of randomness are decoupled. Then, in virtue of 
Eq. (5) we have for the second term 

fdh [ t dt 2 (f(X(t l ),t l )f(X(t 2 ),t 2 )} = ^ ; fdh f dhdiXih) - X{t 2 ))e-^-^' T . 
Jo Jo ^ Jo Jo 

(11) 

Now it is a simple matter to average the delta-function on the rhs of Eq.(ll), by making 
use of its Fourier representation; this yields: 

/°° Ah 
— expHfc(X(ti)-X(f 2 ))) 
-oo 27r 

•exp (—Dk 2 \ti -t 2 \) 



' —CO 

J — c 



' (12) 



yj AllD\t\ t 2 

Here we noted that the process X(t) is Gaussian and we used in the second line the well- 
known rule for averaging Gaussian exponential forms. Inserting Eq.(12) into Eq. (11) and 
integrating leads to: 

<n*)> = *Dt + ((t/r - 1/2) erf [y/t/f] + vWe"^) , (13) 

where erf [x] denotes the error-function, Eq.(7.1.1) in Ref. [14]. Equation (13) is exact 
and is the main result of this section. From Eq. (13) we have in the limit t <C T: 

where the second term represents the result for standard, quenched MdM flows [4]. The 
physics underlying such a behavior is well understood [1,2,4,6,7,9] and it is due to 
dynamically induced correlations in the velocities felt by the particle. 
In the opposite limit, t ^> T, we are lead to 



<y*(*)> = 2Dt + ^ y -t + 0(1) = 2D eff t + 0(1), (15) 
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where 

D <ff = D + ^r)lh> (16) 

i.e. we find diffusion with a renormalized diffusion coefficient. 

Thus for large times t the fluctuations of the environment (flow) contribute additively 
to the diffusion coefficient. Note, however, that now diffusion is anisotropic: Parallel to 
and perpendicular to the layers one has D\\ = D e ff and D± = D, respectively. The 
dependence of D e ff on both D and T is easily visualized: For larger D the perpendicular 
motion is more rapid and hence the velocity field viewed by the particle changes more 
rapidly (and is less efficient); furthermore, a decrease in T leads to a similar effect. On 
the other hand, in the r — > oo limit -D e // diverges, a sign that the motion gets to be 
super diffusive, see Eq. (14). 

Since the disorder changes completely on the time scale T, one may view Eq.(15) 
as arising from renewals of the process in Eq.(14) (with total loss of memory) every r 
units of time, where r is of the order of V. In fact, since the crossover in behavior from 
(2Ar 1 / 2 /3C 2 (7r,D) 1 / 2 )(t/r) 3 / 2 to (A/2C 2 (r,D) 1 / 2 )t occurs at t c = 9nT/lQ w 1.767T, we 
reproduce Eq. (16) exactly if we take r = t c . 



4 Infinitely Long Rouse Chain 



We consider next the opposite case of an infinitely long Rouse chain. For a given thermal 
history, the solution of Eq. (1) for the ^-component of the displacement of the n-th bead 
can be readily found by standard means, see Ref. [1,2], and reads 

Y n (t) = - dr dlP(n-l,t-r)(f(X l (r),r)+r ]y (l,r)), (17) 

S JO J-oo 

where we again, for simplicity, assumed that Y n (0) = 0. Here P(l, r) is the Green's 
function solution of the ID diffusion equation 

d 2 



P(l,r) 



C 



■ cxp 



AKt 



(18) 



A-kKt 

The msd of the n-th bead, averaged both over thermal histories and over random velocity 
fields, takes the form: 

(> \ 1/2 a r-t rt roo poo 

2^Kj + 2TC 2 " Jo dTl I dT2 J dh J dkP{n x 

x P(n -l 2 ,t- r 2 ) e'^- T ^ r 5(X h (r 1 ) - X h (r 2 )) (19) 
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where we used the second moment of the velocity field, given by Eq. (5). 

To proceed further, we need to know 5(Xi 1 (ti) — Xi 2 (r 2 )). This form can be evaluated 
exactly (see Refs. [1, 2, 12, 13] for more details). The result (for t\ < r 2 ) is explicitly: 

5(X h ( Tl ) - X h (r 2 )) = (4^5) V4 M(l 2 - h, n, r 2 )-'l\ (20) 
in which the function M(l, Ti, t 2 ) stands for 

M(/,r 1 ,r 2 ) = (2r 1 ) 1 / 2 + (2r 2 ) 1 / 2 -4(— \ J drP&n + r 2 - 2r) (21) 

Inserting Eqs.(20) and (21) into Eq.(19) and reverting to dimensionless variables, say, by 
setting 

0i = n/t; 6 2 = r 2 /t; 9 = r/t; 

* = &- n imT> * = to- n Km) 1/2 (22) 

we finally arrive at the following expression: 

= 2D (55?) V2 tm + ik (sfe) " 4 (7/49((/r) ' (23) 

with 5(77) being the dimensionless function: 

r02 r°o f oo - z \l (i-0i)-z|/(i-0 2 ) 



/*! /»t/2 /»00 /» OO 

0(77) = / d6 2 dB l e - {ei - d2)r > / d Zl / cfe 2 

Jo 7o J-00 J-00 / 



(1-00(1-02)^ -22,0!,^) 

(24) 

where for 9 2 > 61 the function M(z, 9±, 6 2 ) is 

[61 -z 2 /(9 1 +d 2 -29) 

M{z,0 1 ,e 2 ) = y/W 1 + >JW 2 -2 I d6 , (25) 

JO V ^1 + ^2 — ^ 

In the limit r — > 00 (quenched MdM flows) r] — t/T tends to zero and (7(77) to a constant. 
Hence, for T — > 00 and not-too-small t one has < F n 2 > ~ t 7 / 4 ; this reproduces precisely 
the corresponding result derived in Refs. [1] and [2]. Moreover, we expect this dynamical 
behavior to show up also for finite T, as long as t <C T, i.e. 77 1. 
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Turning now to finite but small T, for f > T, i.e. r\ = t/T ^> f , the leading large-?] 
behavior of g(i]) can be found by integrating in Eq.(24) over 9\ by parts with respect to 
exp(9ir]). We find that gij]) has the form: 



9(r]) = - + 0(-), (26) 
where C is a dimensionless constant, given by 

pi poo poo -(zf+z 2 2 )/(l-e) 

C= d9 d Zl dz 2 = (27) 

Jo J ~°° J -°° {l-9)^M{z 1 - z 2 ,9,9) 

Inserting now Eq.(26) into Eq. (23) we find that the msd of any bead of an infinitely long 
Rouse chain obeys 

««> - 2D {4k) 1,2 (1/2 + f {ikd 1/4 tm + < 28 > 

The significant terms of Eq.(28) can be understood as follows: The first term is just 
the Rouse-result in the absence of flows. The second term is a new feature here, which 
arises due to the time-dynamics of the MdM flows, to be contrasted with the t 7 / 4 -term for 
quenched MdM flows. This new term appears, of course, even when the time-correlation 
in Eq. (5) is of <5-form, i.e. for r — > 0. 

Therefore, we have for monomers attached to long polymer chains that at long times 
< Y£{t) >~ t 3 / 4 , i.e. that the exponent of the anomalous diffusion is larger than in the 
force-free case. One may contrast this behavior to that of a single particle, where at 
long times time-dependent velocity fields only change the diffusion constant, but leave 
the exponent of t (unity) unchanged. Hence the influence of the velocity fields on the 
monomer's motion differs from that of the thermal noise. The difference can be traced to 
the form of the disorder; the solvent's role is included through MdM 5(Xi 1 {ti) — X^(t 2 ))- 
fields, which depend on the X-components of the positions of the beads, while thermal 
noise depends on Si u i 2 . This means that while different monomers always experience 
different forces from the heat bath, distinct monomers with the same X-components 
(there are many such monomers, since the polymer's projection on the X-axis corresponds 
to a simple random-walk) experience identical velocity fields, a fact which qualitatively 
enlarges < Y^(t) >. 
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5 Finite Rouse Chains 



For finite Rouse chains there appears yet another timescale, namely the so-called Rouse 
time tn = (N 2 /(ir 2 K) [11-13], t R being the largest internal relaxation time of the struc- 
ture. To display its role we start by solving Eq. (1) for finite N. The solution is found 
most readily in terms of normal coordinates [12], introduced, say for Y n (t), through: 

oo 

Y n (t)= cos(pnn/N)Y(p,t), (29) 

p=— oo 

and similarly for X n (t) and Z n (t). Inserting these forms into the equation of motion, 
Eq. (1), and solving yields straightforwardly for the X-component [1] 

n ft pN 

X(p,t) = — I dr dnr] x {n,r)cos(pnn/N)e- p2{t - T),tR (30) 
Jo Jo 



and similarly for the Z-component. For the ^-component one has instead 
1 f l f N 

Y M = TF?I dT dn Mn, r) + f(X n (r),r))cos( I mn/N)e-^ t - T ^. 
iV< 5 Jo Jo 



(31) 



The behavior of the beads depends on the ratio t/t R . For short times t <^t R the sums over 
p for X n {t), Y n (t) and Z n (t), e.g. Eq. (29), can be converted into integrals; exemplarily, 
inserting Eq. (31) into Eq. (29) and performing the integration leads to Eq. (17). Thus 
for t <C t R each bead behaves as if it were part of an infinite chain: Hence the results 
from the previous section apply. In the opposite regime, for t ^> tR, only the zeroth mode 
p = contributes significantly to Eq. (31), and the motion of each bead follows closely 
that of the center of mass (CM) of the chain. In the case of the X- and Z-components, 
this is a diffusive behavior with the renormalized diffusion constant D — > D/N [11-13]. 
For the F-component, however, we have for t~^> t R approximately: 

D 



(Y 2 (t))^(Y 2 (0,t)) = 2-t 



A 

+ 



2N 2 T(- 



r-t ft i-N fN 

dn dr 2 dm dn 2 5(X ni (n) - X n2 (r 2 )) e -^- T ^ r . (32) 
Jo Jo Jo Jo 



We remark that the terms on the rhs of Eq. (32) describe the dynamics of the ^-component 
of the chain's CM for all t. The evaluation of Eq. (32) is rendered complex for finite N due 
to the appearance of the two-time delta function in the integral. It is, however, possible 
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to proceed along the lines of Ref. [2], but this is outside the scope of the present article. 
Since the analysis simplifies for r — > oo and for r — > 0, we shall consider only these two 
cases in the following. 

For r — > oo the flows get to be time- independent; then the approach follows the 
derivation given in [1, 2] for quenched MdM flows. There it was found that the leading 
behavior of the bead's motion obeys (Eq. (36a) of Ref. [1], here in our notation): 

(W)) ~ r^V^ 372 ' (33) 

that is, the CM moves practically as an individual particle, i.e. Eq. (14), with D being 
however replaced by D/N. 

For r — > 0, we have for the second term (call it /) in Eq. (32) 



pt pN pN 

dr dm dn 2 S(X ni (r)-X n2 (r)). (34) 
Jo Jo Jo 



NX 1 

For t^>t R most of the r in Eq. (34) also obey r > so that we may view the Gaussian 

X-process to be stationary, i.e. X ni (r) — X m (r) = -D£l n i — n 2 \/K, independent of r. 
Hence 

S(X ni (r) - X n2 (r)) = / WM-^W) 

Jo 21T 

°° dk k 2 DC\ni-n 2 \/2K 



o 2vr e 



K 



2D7r£|ni — n 2 \ 

Inserting now Eq. (35) into Eq. (34) and performing the integrations leads to 



(35) 



A IK r , f N , f N , , ,_ 1/2 8A / K t 

dr I dn\ / dn 2 \ri\ — n 2 \ ' =-—\l————= (36) 



NX V 2DttC Jo Jo Jo ' ' 3C 2 V yfN 



Using this expression in Eq. (32) we see that (Y%(t)) = D^fft, where now -Dg/j is 



D* - D + 4A K 1 
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Again the time dependent MdM flows are seen to produce a "correction" to the diffusion 
constant at long times; this correction is of more importance the longer the chains. For 
infinite chains vanishes, a sign of the appearance for N — > oo of the sublinear, t 3 / 4 - 
behavior of Eq. (28). The result can be understood in terms of the renewal of the process 
of Eq. (28): Thus for t 3> t R 3> T, assuming a loss of memory after each t R units of time, 
(Y£{t)) can be expressed as being around t/t R times C (n( 2 )(K / (W^Cf/H^ . This 
leads to a correction of the diffusion constant of 

j_ca fj^v 1 ' 3/4 = nr i 

2t R vrC 2 \4DH() " 87rV4 3C 2 y 2D7TC ViV 1 ' 

which is 3C/(87r 1 / 4 ) times the correction in Eq. (37). 

A similar argument can be put forth when t ^> F ^> t R . In this case, (Y%(t)) can be 
viewed as arising from renewals (every T units of time) of the process given by Eq. (33). 
This leads to yet another effective diffusion constant D^ff, which besides D/N has an 
additional term proportional to 

a rw_ 

This reproduces (after replacing D/N by D) the result of the single particle case, Eq. (16). 

We summarize the findings for the msd of a bead on the finite Rouse chain. When 
t ^> max(t R ,T), there is normal but anisotropic diffusion of the bead under observation, 
described by the diffusion constant D± = D in the directions perpendicular to the flow, 
and an effective diffusion constant D\\> D parallel to the flow. The value of D\\ depends 

on whether t R < T or t R > T: For t R <^ T, D\\ = D^ff ( see E Q- ( 39 )) and for *fl > r > 
D\\ = D^jp Eq. (37). In the short time regime, t min(^,r), we always have Y%(t) ~ 
t 7 / 4 . However, the intermediate regime again depends on the magnitudes of t R and T. If 
t R > T we observe (Y r ^(t)) ~ t 3/2 (Eq. (33)) in the regime t R < t < T. But if on t he other 
hand T 3> t R) then for times t obeying T t t R , the msd behaves as (Y%(t)) ~ t 3//4 . 
Our model gives therefore rise to a very rich dynamical behavior. 

6 Conclusion 

In this work we introduced a new variant of the Matheron-de Marsily (MdM) model of 
layered random media, in which the velocity field is also allowed to change with time. 
We solved exactly the equations describing the dynamics of a single particle and of beads 
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(39) 



belonging to infinite Rouse chains in time-dependent MdM flows. For short times and 
slowly changing media, we recovered the results of the standard quenched MdM model. 
For large times we found normal diffusion in the cases of a single particle and of a fi- 
nite Rouse chain, and we computed the corrections to the diffusion constant. For the 
infinite Rouse chain we found a surprisingly fast increase with time of the mean square 
displacement, a result also valid for finite Rouse chains in intermediate time regimes. 
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